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ABSTRACT 

We report the discovery of 8.5cr high-frequency quasi-periodic oscillations (HFQPOs) at 66 Hz in the 
RXTE data of the black hole candidate ICR J17091-3624, a system whose X-ray properties are very 
similar to those of microquasar GRS 1915+105. The centroid frequency of the strongest peak is —66 
Hz, its quality factor above 5 and its rms is between 4 and 10%. We found a possible additional peak 
at 164 Hz when selecting a subset of data; however, at 4.5c level we consider this detection marginal. 
These QPOs have hard spectrum and are stronger in observations performed between September and 
October 2011, during which ICR J17091-3624 displayed for the first time light curves which resemble 
those of the 7 variability class in GRS 1915+105. We find that the 66 Hz QPO is also present in 
previous observations (4.5er), but only when averaging —235 ksec of relatively high count rate data. 
The fact that the HFQPOs frequency in IGR J17091-3624 matches surprisingly well that seen in 
GRS 1915+105 raises questions on the mass scaling of QPOs frequency in these two systems. We 
discuss some possible interpretations, however, they all strongly depend on the distance and mass of 
IGR J17091-3624, both completely unconstrained today. 

Subject headings: X-rays: binaries — binaries: close — stars: individual (IGR J17091-3624, GRS 
1915+105) — Black hole physics 



1. INTRODUCTION 

One of the strongest motivations for studying Low- 
mass X-ray binaries (LMXBs) has been the aim to use 
these systems as probes of fundamental physics. It 
has been argued that the study of high frequency (40- 
450 Hz) quasi-periodic oscillations (HFQPOs) found in 
LMXBs containing a stellar mass black hole (BH) is 
one of the most promising tools to do so. The HFQ- 
POs are the fastest phenomenon observed in these sys- 
tems and their high frequency suggests that they are 
produced in the innermost region of the accretion flow 
near the BH event horizon. Once the processes in- 
volved in their formation are firmly understood, they 
will give a wealth of information on BH mass and 
spin as w ell as the structure of strongly curved space- 
time (e .g..lNowak"et^ll^9^ r Wagonedll99llStella et alJ 



1999; 



2001 



iWagoner et alJ l2001t lAbramowicz fc Kluzniakl 
Katdl2001t IRezzoUa et al.ll2003l ~ 

The frequency of the HFQPOs usually occur at spe- 
cific values, different in each source. In a few cases, pairs 
of QPOs have been detected. The highest frequency 
appears to scale inversely proportional to the to BH 
mass (e.g., iMcClintockl ^ Rcmillard 2006; iBelloni et all 
20063- HFQPOs are weak, transient, and energy depen- 
dent, and are mostly detected only at high count rates. 
In some cases they are found just above detection lev- 
els, so there is considerable uncertainty about their ex- 
act properties. However, typically the quality factor Q 
ranges between less than 2 and a few tens, and ampli- 
tudes are between 0.5 and 5% rms (depending on the 
source and the energy range used) . 
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Twi n QPOs at a stable 300 Hz and 450 Hz (|Strohmaverl 
I2001af) have been observed in t he BH GRO J1655-4 0; at 
240 and 160 Hz in H1743-322 (|Homan et al.l[2005T ) and 
at 188 and 268 Hz in XTE J1550-56 (note however, that 
in this case it is n ot clear how stable th e frequency of 
these QPOs is, se e Rcmill ard et al.l Il999t iHoman et alJ 
120011 ; iMiller et al.ll2001h . Single QPOs or peaked noise 
(i.e. when Q <2) features at a rather stable frequency 
have been detected in XTE J1650-500 (-250 Hz, see, e.g, 
Homan et al.ll2003f). 4U 1630-47 (variable frequency, see 



Klein- Wolt et al J 12004) a nd XTE J1859+226 (-180 Hz, 
see, e.g.. ICui et al.ll2000l ). 

The remaining BH which shows HFQPOs is the pe- 
culiar system GRS 1915+105, which is a 33.5 days or- 
bital perwcHaixu^ Mq 
BH (jGreiner et alj[200ll; lHarlaftis fc Greinerl [200l At 
a distance of — 12.5 kpc ((Mirabel fc Rodriguezl ll994) . it 



is ver y often at Eddingt on or super-Eddington luminosity 
(e.g./ lDone eTall I2004D . Until recently, GRS 1915+105 
has been unique in that its X-ray light curves ex- 
hibit more than a dozen different patterns of variabil- 
ity usually called "classes" (which are referred to with 
Greek letters), most of which are high am plitude and 
highly-structured (e.g., IBelloni et alJ |2000T ). Most of 
this structured variability occurs on timescales of sec- 
onds or longer, and it is thought to be due to limit cy- 
cles of accret i on an d ejection in an unstable disk (e.g., 
I Belloni et ail 119971: iMirabel et alJ 119981 : iTagger et all 
12004 iNeilsen et alJl20lif T 

GRS 1915+105 has HFQPOs at the r ather low fre- 
quencies of -40 Hz (jStrohmaverl |2001H ) and -67 Hz 
(|Morgan e t al. 1997) compared with the rest of the BH 
systems; however, as in the case of GRO J1655-40 and 
H1743-322, the QPO frequencies remain almost constant 
as the X-ray flux changes (any change in the frequency 
of the —40 Hz oscillation during the sequence of obser- 
vations was <1%; the frequency of the 67 Hz oscilla- 
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Fig. 1. — Upper panel: Crab-normalized 2-20 keV intensity of the full outburst of IGR J17091-3624 sampled with RXTE. Lower panel: 
0.01-10 Hz frac tional rms amplitud e (2-60 keV) range. Points are per observation. The grey area marks the period of fast 7 class of 
variability (e.g., Bclloni ct al. 2000). The average power spectra of all observations in this region is shown in Figure[2] 



tion varies by only a few percent as the X-ray flux varies 
by a factor of s everal) , suggesting; a common origin be- 
tween systems. iBelloni et al.l ()2001l ) reported on the de- 
tectio n of QPOs at ~27 Hz; given the weakness of the 
signal Be lloni et ail (|200lD could not determine whether 
the ~27 Hz was related to the -67 Hz QPO (at the 
time of IBelloni et al.l 120011 the -40 Hz QPO was not 
yet reported); however, they were able to conclude that 
the presence/absence of HFQPOs is intimately linked 
to the slower oscillations and variations that happen on 
timescales longer than seconds. An additional peak at 
~170 Hz was also found after ap plying a very specifi c 
class-count rate- hardness selection ([BeHoni et al.l l2006b); 
however, this barely significant peak is very broad, with 
a Q value or —2. 

If the frequency of the HFQPOs depends on the mass 
and/or spin of the BH, one would expect (as observations 
appear to suggest), that different systems with differ- 
ent parameters should show QPOs at different frequen- 
cies. The recent discovery of a BH (IGR JI7091-3624) 
which shows very similar high amplitude and highly- 
structured variability to that seen in GRS 1915+105 (see 
lAltamirano et aLll2011cL and references therein) opened 
a new window of opportunities to understand the physi- 
cal mechanism that produces the highly structured X-ray 
variability and it relation with much faster (sub-second) 
variability. 

IGR J17091-3624 has been seen active in 1994, 



1996 2001, 2003, 20 7 and 2011 (iRevnivtsev et all 
20031: lin't Zand et all 120031: ICapitanio et al.l 120061: 



Kuulkers et al.ll2003t ICapitanio et alll2009l: iKrimm et all 
2011J). Only hi the recent 2011 outburst IGR J17091- 



3624 started to show t he high amplitude and h ighly- 
structured variab ility (jAltamirano et all l2011dH lbl|a: 
iPahari et al.l 121311 typical of GRS 1915 +105. In the 
contex t of the variability class es defined b y Bcllon i et all 
j2003) for GRS 1915+105, lAltamirano eTail ^OlUh 
found that IGR J17091-3624 shows the v, p, a, A, 
P and fi classes as well as quiet periods which re- 
semble the x class, all occurring at 2-60 keV count 
rate levels which can be 10 -50 times lower t h an ob - 
served in GRS 1915+105. lAltamirano et all (|2011cD 
also found that the difference in flux combined with 
the circumstance that currently neither the distance 
to IGR J17091-3624 nor the mass of its compact 
object are known lead to the conclusion that either 
all models requiring near Eddington luminosities for 
GRS 1915+105-like variability fail, or IGR J17091-3624 
lies at a distance well in excess of 20 kpc or, it harbors 
one of the least massive BHs known (< 3M Q ). 

The similarities and differences between IGR J17091- 
3624 and GRS 1915+105 raises the question of at 
what frequency should we observe the HFQPOs in 
IGR J17091-3624. Assuming 1/M scaling, and the pos- 
sibility that IGR J 1709 1-3624 has a mass which is a fac- 
tor of -5 smaller than GRS 1915+105, then the HFQ- 
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Fig. 2. — Upper panel: Averaged power spectra of all observa- 
tions of IGR J17091-3624 marked by the grey area in Fig. [T] (for 
a total of ~74 ksec of data). Lower panel: Average power spectra 
of a single observation (ObsID: 80701-01-28-01) of GRS 1915+105. 
Both plots show the typical 7 variability class power spectra, with 
a strong broad bump at low frequencies. Insets show a zoom to 
the ~ 65 - 67 Hz QPO. 

POs in IGR J17091-3624 should be at a frequency sig- 
nificantly higher than the ~40 Hz and ~67 Hz seen in 
GRS 1915+105. Finding similar HFQPO frequency be- 
tween sources would imply that the frequency does not 
depend on the mass or that IGR J17091-3624 lies at a 
distance well in excess of 20 kpc. 

Triggered by this question, we searched all the 
Rossi X-ray Timing Explorer (RXTE) available data 
of IGR J17091-3624 for the presence of QPOs at 
frequencies higher than the ~7-10 Hz reported by 
lAltamirano et all (|2011cf ). In this Letter we report the 
discovery of a highly significant QPO at ~65 Hz, and a 
marginally detection of a QPO at ^165 Hz and briefly 
discuss the implication of our detections. 

2. OBSERVATIONS AND DATA ANALYSIS 

IGR J17091-3624 was observed with the P roportional 
Counter Array (PCA; iJahoda et al.1 [2006) on-board 
RXTE almost daily since the outburst began in Febru- 
ary 2011 until October 15th, 2011, when RXTE could 
not point to IGR J17091-3624 any more due to visibility 
constraints. A total of 228 observations, covering ~ 280 
days, sample the 2011 part of the current outburst (At 
the moment of submission it is not known whether the 
source is still active or not). In this work we do not use 
data taken between MJD 55600 and 55614 (14 observa- 
tions), given that these observations were affected by a 
bright source in the field of view (e.g., lAltamirano et all 
I2011cf ). 

We use the 16-s time-resolution Standard 2 mode data 



to calculate the Crab-normalized ave rage per observa- 
tion 2-20 keV intensity (see, e.g., lAltamirano et all 
l2008l and references therein for details on the method). 
Power spectra and light curves were produced from the 
PCA using standard te chniques (e.g.. lBelloni et al1l200Ct 
lAltamirano et aLll20 08). For the Fourier timing analysis 
we used data in the 2-25 keV range (absolute channels 
5-60) from the Good Xenon mode. 

To fit the high-frequency (> 10 Hz) part of power 
spectra, we used one or two Lorentzian functions plus 
a power-law when needed. We give the frequency of the 
Lorentzians in terms of characteristic frequency v max = 
y^n + (FWHM/2) 2 = i/ ov /l + 1/4Q 2 (jBeUoni et al.l 
2002). For the quality factor Q we use the standard 
definition Q = v /FWHM. FWHM is the full width 
at half maximum and v$ the centroid frequency of the 
Lorentzian. The strength of the QPOs are given in terms 
of fractional rms amplitude. 

3. RESULTS 

3.1. Outburst evolution 

In the upper and lower panel of Fig. [TJ we show the 
average 2-20 keV intensity for each observation and the 
0.01-10 Hz rms amplitude, respectively. During the 265 
days period, IGR J 1709 1-3624 has been at an aver- 
age intensity between 20 and 65 mCrab. (Note that 
during single flares or heartbeats the source was much 
brighter). The rms amplitude also varied significantly, 
where IGR J17091-3624 underwent periods of very low 
(~10%) and very high (>60%) amplitude. Figure Q] 
shows that there is no clear correlation between intensity 
and rms amplitude. A similar study of GRS 1915+105 
shows the same (this will be reported elsewhere). 

3.2. High frequency QPOs 

At first we searched each observation power spectrum 
for indications of QPOs at higher frequency than 10 Hz 
and at a significance larger than 3.5<r. We found a clear 
detection (4.7<r, single trial) of a QPO at v — 64.8 ± 
1.7 Hz, Q = 4.8 ± 1.7 and rms amplitude of 10 ± 1% 
(2-25 keV, ObsID: 96420-01-36-03, MJD 55865.5). We 
also found that in the period MJD 55830-55880 many 
power spectra revealed power excesses in the 60-70 Hz 
range. To further investigate this, we added all data 
in this period , which corresponds to the low-variability 
period marked with a grey area in Fig. [TJ In the upper 
panel of Fig. [2] we show that the resulting power spectrum 
is clearly dominated by a strong feature at ~0.18 Hz, but 
also shows an 8.5cr QPO at 66.5 ± 0.5 Hz (Q = 7.8 ± 1.5; 
4.9 ± 0.2% rms amplitude). 

Different sub-selections on time and intensity of these 
data lead to significant changes in Q, amplitude and sig- 
nificance, but not in frequency. This suggests that that 
the 66 Hz QPO is not always present, similar to what oc- 
curs in GRS 1915+105. In several of our sub-selections 
we saw an excess of power in the 150-200 Hz range. In 
the left panel of Fig. [3] we show the result of averaging all 
observations after MJD J55841. The QPO at - 66 Hz 
(9a significant, Q = 6 ± 1, 6.4 ± 0.4% rms amplitude) 
is accompanied with a a QPO at 164 ± 10 Hz (~4.5cr, 
Q = 2.6 ±1.3, 6.6 ±0.8% rms amplitude); In the inset we 
show the rms amplitude vs. energy of the 66 Hz QPO; 
its amplitude increases with energy from ~5% at 3 keV, 
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Fig. 3. — Left panel: Average power spectrum of all observations after MJD 55841 (for a total of 48 ksec). The 66 ± 0.5 Hz QPO is 9a 
significant, while the ~ 164 ± 10 Hz is 4.5<r. Lower inset shows the rms amplitude vs. energy for the 66 Hz QPO. Bight panel: Average 
power spectrum of all observations performed before MJD 55830, when the average intensity of IGR J17091-324 was higher than 46 mCrab 
(corresponding to ~100 cts s _1 PCU -1 ; we used a total of ~235 ksec). Note the difference the y axis scales. The broad (Q=l) QPO at 
70 ± 6 is 4.5<r significant. In both panels we used data in the 2-25 keV range. Insets on the top: light curves marking with grey area the 
data used to produce each power spectrum. 

to -14% at 13 keV. The 67 Hz QPO in GRS 1915+105 
also has a hard spectrum , although less extreme than this 
one iMorgan et alj (|1997D . The low S/N did not allow a 
similar analysis for the 164 Hz QPO. 

We also explored whether the 65 Hz QPO is also 
present before MJD 55830 which is not detected due 
to the low-statistics. Averaging all data in this period 
led to a broad feature at —60-80 Hz whose parameters 
were unconstrained. When we selected only high count 
rates (> 110 cts/sec/PCU), we detected a 4.5cr broad 
(Q = 1.2 ± 0.5) feature at 70 ± 6 Hz (4.4 ± 0.4% rms 
amplitude). We show our best fit n the right panel of 
Fig. El 

We did not detect any QPO at 40 Hz nor 27 Hz as 
seen in GRS 1915+105. This is not surprising, as those 
QPOs are can be weaker and com plex data selection 
might be needed t o find them (e.g.. IStrohmaverll2001bl : 
IBelloni et al1[200l 

3.3. The low-variability period MJD 55830-55880 

The period MJD 55830-55880 (where the 66 Hz QPO 
is strongest, see Section [3~2]) corresponds to an interval 
of relative low (10-30%) rms amplitude. The 1-sec light 
curves during this period are relatively flat, apart from 
quasi-periodic oscillations with a typical time scale of 10 s 
and/or the presence of sharp dips with a typical duration 
of a few seconds. In the upper panel of Fig. [4] we show a 
representative lightcuve during this period. 

When compare d with the GRS 19 15+105 variability 
classes defined bv IBelloni et~aT1 ()2000t ). the IGR J17091- 
3624 light curves during MJD 55830-55880 are reminis- 
cent of the 7 class. As a comparison, in the lower panel 
of Fig. [4] we plot a representative 7 class data-segment 
for GRS 1915+ 105 (for the full light curve, see Figs. 
2f and 15b,d in IBelloni etall [2000) . A comparison be- 
tween power spectra between the IGR J 1709 1-3624 and 
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Fig. 4. — Upper and lower panel show 1 sec bin light 
curves of IGR J17091-3624 (subset of ObsID 96420-01-37-04) and 
GRS 1915+105 (s ubset of observation 20402-01-39-00, see also 
Bclloni et al. 2000), respectively, during the 7 variability class. 

GRS 1915+105 is shown in Fig. [2 

4. DISCUSSION 

We discovered the presence of a HFQPO in the RXTE 
data of the IGR J17091-3624, a system whose X-ray 
properties are very similar to those of GRS 1915+105. 
The centroid frequency of the strongest peak is —66 Hz, 
its quality factor above 5 and its rms is between 4 and 
10%. We found a possible additional peak at 164 Hz 
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when selecting a subset of data; however, we consider 
this detection marginal. 

The HFQPOs are detected during a period where 
IGR J17091-3624 i s in the 7 variabili ty class 
(jBelloni et all l2000f ) . lAltamirano et al.1 l2011d reported 
observations until MJD 55800 where so far IGR J17091- 
3624 was seen in at least the v, p, a, A, /3, (i and x 
classes. The fact that now we also identified the 7 class 
probably means that IGR J17091-3624 is still evolving 
between classes. 

The frequency of the 66 Hz QPO in IGR J17091-3624 
is surprisingly similar to that of the main HFQP Q in 
GRS 1915+105 (65-69 Hz, see IMorgan eTaTl \l99% . In 
addition, the frequency of the 4.5ct significant broad fea- 
ture at 164 ± 10 Hz in IGR J17091-3624 also matches 
that of a similar broad featu r e foun d at 170 Hz in 
GRS 1915+105 (|Belloni et all l2006al) . The circum- 
stance that currently neither the distance to IGR J17091- 
362 4 nor the mass o f its c ompact object are known 
led lAltamirano et all (|2011d) to the conclusion that ei- 
ther all models requiring near Eddington luminosities for 
GRS 1915+105-like variability fail, or IGR J17091-3624 
lies at a distance well in excess of 20 kpc or, it harbors one 
of the least massive BHs known (< 3M ). The surpris- 
ingly match between the ^70 Hz and ^170 Hz frequency 
therefore raises more questions than answers. 

The frequencies of dynamical motions n ear compact 
objec ts are thought to scale with mass (e.g.. Ivan der Klisl 
2006). The slower highly-structured low- frequency 
variability seen in GRS 1915+105 as compared with 
IGR J17091- 3624 supported the idea that the later was 
less massive ()Altamirano et aLll2011d ). However, the fact 
that the HFQPOs occur at the same frequency could be 
at odds, meaning that either the frequency of the HFQ- 
POs or the highest-frequency of the 'heartbeats' in these 
two sources scale with mass, or that none of them do, or 
that they both do thanks to a difference in BH spin. 

If the 67 Hz HFQPOs frequency scales with mass (as- 
suming same mass, same spin), then either IGR J17091- 
3624 is further than 20 kpc or it is close, in which 
case all models requiring near Eddi ngton luminosities 
for GR S 1915+105-like variability fail (jAltamirano et ail 
l2011d) . Under the assumption that the centroid frequen- 
cies represent that of a Keplerian circular motion at the 
inner-most stable orbit (ISCO), then at first approxima- 
tion (e.g., iKluzniak et al.l [ i~9 90) no solution is possible 
for the 67 Hz QPOs independently of the spin, unless 
either the HFQPOs are not produced at the ISCO or 
the mass of both BHs is significantly higher than the 
~14M estimated for GRS 1915+105. The 160-170 Hz 
broad QPOs found in both in IGR J17091-3624 (this 
Letter) or GRS 1915+105 would be consistent with a 
14-2OM , but only if the spin parameter (o/Mbh) in 
these systems is significantly lower than 1 (i.e., at odds 
with the reported a/M B H ~ 1 for GRS 1915+105, e.g., 
IMcClintock et ah 2006; Bl um et alj[2009h . 



In th e framework of th e Relativistic Precession Model 
(RPM. [Stella et aI1ll999D . and relaxing the assumption 
on equal spin, we can identify the higher frequency (160- 
170 Hz) as representing a Keplerian motion at a certain 
radius and, the lower frequency (66 Hz) as the relativis- 
tic periastron precession at the same radius (notice that 
this radius must be larger than the ISCO, where the two 
frequencies are the same). To test this identification, 
we searched for the combinations between spin parame- 
ter and BH mass which allow a solution. We find solu- 
tions for masses ranging between 7 M (for cl/Mbh = 0) 
and 19 M (for o/Mbh = 1); the latter consistent 
with the 14.0 ± 4.4 M w and a/M B H — 1 reported for 
GRS 1915+105 dGreiner et alJl200lt lHarlaftis fc Greinerl 
[200l IMcClintock et al.l 120061: fBlum et al.1 12009ft . This 
identification implies that IGR J 1709 1-3624 could be 
considerably less massive than GRS 1915+105 (and coin- 
cidentally have a much lower spin than GRS 1915+105), 
explaining in this way why the highest frequency of the 
'heartbeats' is significant ly higher in IGR J17091 -3624 
than in GRS 1915+105 (jAltamirano et al.ll2011d) . We 
note that in the framework of the Relativistic Precession 
Model the 41 Hz frequency would be unexplained. Solu- 
tions to the 41,67 Hz pair for GRS 1915+105 are possible, 
but only for masses larger than 25M and considering no 
spin (i.e o/Mbh = 0). 

If the HFQPOs frequency does not scale with mass, 
then it is intriguing why their HFQPOs frequencies are 
the same in these two sources. One possibility is that the 
67 Hz HFQPOs are different from the HFQPOs seen in 
other sources (e.g., Ivan der Klis 200j|), and that their 
frequency is instead related with the physical process 
that produces the highly-structured low-frequency vari- 
ability (only seen) in these two sources. If true, this 
could have major implications on works that use the 
~ 67 Hz HFQPO in GRS 1915+105 as a standard to 
estim ate the mass of galactic a nd/or super massive BHs 
(e.g. iMiddleton fe Dondl20Tfl . 

Our results add another HFQPO to the small num- 
ber already known. Its properties are not different from 
the others, in particular its elusiveness. There are thou- 
sands of observations of more than 20 black-hole tran- 
sients in the RXTE archive, but only very few signifi- 
cant peaks have been found. As RXTE will cease op- 
erations in 2012 January, this sample is not going to 
be increased until the Indian astronomical mission AS- 
TROSAT is launched, while a major step forward will 
only come with t he next generation of X-ray instruments 
(e.g. LOFT, see lFeroci et~aTll2011h . 
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